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A deep X-ray observation of supernova remnant 
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ABSTRACT 

In this paper, we present the analysis of a deep (99.6 ks) observation of G304. 6+0.1 
with the X-ray Imaging Spectrometer on board Suzaku satellite. The X-ray spectral 
data are well-fitted with a plasma model consisting of a thermal component in colli- 
sional ionization equilibrium and a non-thermal component. The thermal emission is 
well fitted with VMEKAL model with an electron temperature of kT e ~ 0.75 keV, a 
high absorbing column density of JVh ~ 3.9 x 10 22 cm~ 2 and near/lower solar abun- 
dances which indicate that the X-ray emitting plasma of G304. 6+0.1 is dominated by 
swept-up ambient medium. The non-thermal component is well fitted with a power- 
law model with photon index of T ~ 1.4. We found a relatively high electron density 
n c ~ 2.3/~ 1/2 cm" 3 , age t ~ 1.4 X 10 4 / 1/2 yr, and X-ray emitting mass M x ~ 380/ 1/2 
M© at an adopted distance of d=10 kpc. Using the morphological and spectral X-ray 
data, we confirm that the remnant is a new member of mixed-morphology supernova 
remnants. 
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1 INTRODUCTION 

Recent supernova remnant (SNR) observations in X-ray 
band led to discovery of a considerable number of a new 
morphological class called mixed-morphology (MM) SNRs 
that have a centrally bright in the X-ray band characterized 
by a thermal spectrum with little or no limb brightening 
and a shell morphology in radio ( Rho fc Petre|1998 I . Their 
characteristics have been explained with two basic scenar- 



(19701. Caswell et al. (19751 gave a lower limit of 9.7 kpc 



to the distance of the remnant by measuring 21 cm absorp- 
tion profiles with the Parkes hydrogen line interferometer. 
From the radio observations, Milne et al. (1985) suggested 
that its structure has been resulted from the shock front 
collision with its dense environment, Whitcoak fc Green] 
(1996) have shown that it had a irregular shell morphol- 



ogy with a small size. Frail et al. (J1996J) have detected 



ios; evaporating clouds (Cowie & Mckee 1977 White & Long 



1991) and thermal conduction models (Cox et al. 1999 Shel 



ton et al.|1 999 Shclt on, Kuntz fc Petre|2004| ). In addition to 
these basic ones there are new scenarios that try to explain 
the mechanism taking place, but MM SNRs are still not un- 
derstood completely. Deep X-ray observations will help un- 
derstanding the mechanism taking place in these MM SNRs. 
G304.6+0.1 (Kes 17), a new MM member (|Combi et 



al.||2010| ), was first detected at 408 MHz and 5 GHz with a 
Parkes and a Molonglo radio telescopes by Shaver & Goss 
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1720 MHz OH maser e mission toward G304. 6+0.1 by us- 
ing a Parkes Telescope. |Reach et al.| ( |2006[ ), using mid- 
infrared data from the Spitzer Space Telescope at 3.6, 4.5, 
5.8 and 8 /im showed that the remnant has a shell structure 
in infrared and detected a pair of thin filament at about 
RA(2000) = 13 h 05 m 46 s .2, Dec. (2000) = -62°38'33". The 
authors claimed that most of the infrared emission com- 
ing from the shell could be resulted from shocked molecular 
material. Hewitt et al. (20091 have obtained sensitive low- 



resolution spectroscopy from 5 to 35 /im for this remnant 
with Spitzer IRS and showed evidence of shocked molec- 
ular gas. The first detailed multi-wavelength (radio, mid- 
infrared, and X-ray) study of the G304. 6+0.1 has been con- 
ducted by |Combi et al.| ( |2010) . By using XMM-Newton data 
(26.8 ks) they have found that the X-ray emission was dom- 
inated by thermal radiation (PSHOCK) plus a small frac- 
tion of non-thermal radiation (power-law) in north, center 
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and south regions. They have found a high column den- 
sity iV H in the range of (2.5 - 3.5) xlO 22 cm" 2 , kT c ~ 0.75 
keV and plasma was in non-equilibrium ionization state with 
r ~ 2.1 x 10 12 cm~ 3 s. In addition, they calculated elec- 
tron density n c to be ~0.99, ~0.89 and ~2.26 cm" 3 for 
north, center and south regions, respectively. Their multi- 
wavelength study has suggested that this remnant was a 
middle-aged SNR, and belonged to the MM SNRs class. |Lee| 



et al. (2011 1 showed that the IR emission was concentrated 



on the western and southern shells by using data obtained 
with AKARI and Spitzer satellites, the western shell was 
bright in the mid and far-IR continuum and the near-IR H2 
line emission, whereas the southern shell was visible only 
in the mid-IR continuum emission. The authors suggested 
that far-IR continuum and near-IR H2 line emission emanat- 
ing from the western shell was resulted from the interaction 



with nearby molecular cloud. Wu et al. (20111 reported the 



discovery of GeV emission at the position of this remnant 
by using the data from the Large Area Telescope on board 
the Fermi Gamma-ray Space Telescope. They observed the 
remnant in the 1—20 GeV range with a significance of ~12<r, 
and suggested that Kes 17 was a candidate acceleration site 
for cosmic-rays, and reported the detection of a number of 
7-ray sources in its vicinity. 

Suzaku ( Mitsuda et al.||2007 l, the newest Japanese X- 
ray astronomy satellite that has a good energy resolution, 
high efficiency, and low background level for diffuse sources, 
has observed G304. 6+0.1 with 99.6 ks exposure time. Our 
aim is to study the remnant whose morphology implies a 
similarity to the recently proposed group of MM SNRs by 
using the high quality spectral data. 

The rest of this paper is structured as follows. Obser- 
vation and data reduction are described in Section 2, image 
and spectral analysis are reported in Section 3. We discuss 
results of the spectral analysis and finally, investigate the X- 
ray characteristics of G304. 6+0.1 from radial spectral varia- 
tions and compared the results with the two basic scenarios 
explaining MM SNRs in Section 4. 



2 OBSERVATION AND DATA REDUCTION 

Suzaku observed G304. 6+0.1 with the X-ray imaging spec- 
trometer (XIS: |Koyama et al.| ( |2007[ )) on 2010 September 
03, for 99.6 ks (Obs ID:505074010). Suzaku has three work- 
ing XIS CCDs, two of the cameras (XIS0 and XIS3) have 
front-illuminated (FI) CCDs, and the remaining one (XIS1) 
has a back-illuminated (BI) CCD. Each CCD camera has a 
single CCD chip with an array of 1024 x 1024 pixels, and 
covers an area of 17.8x17.8 arcmin 2 . The XIS was operated 
with normal clocking mode, in 5 x 5 and 3x3 editing modes. 
We used the HEADAS 6.5 sofware package for data reduc- 
tion. The XIS data have been screened with xselect using 
the standard criterion ( jKoyama et al.|2007[ ). Spectral analy- 
sis and model fitting have been performed with xspec vll.3 



( Arnaud|[l996| ). The redistribution matrix files (RMFs) of 
the XIS has been produced by xisrmfgen and auxiliary re- 
sponce files (ARFs) by xissimarfgen (Ishisaki et al.|2007 1. 



3 ANALYSIS 

3.1 Image analysis 

Fig. 1 shows XIS0 image of G304. 6+0.1 in 0.3-10 keV full 
energy band. The outer most ellipse represents the over all 
region (we call region 1) that X-ray emission originates. In- 
side the region 1, the emission looks much concentrated in 
two parts, inner ellipse (we call region 2) and inner circle 
(we call region 3). The dashed circle and dashed square show 
background region and the XIS field of view (FOV) , respec- 
tively. Fig. 2 shows the XIS0 image in 0.3—10 keV band 
overlaid with contours of the 843 MHz radio image from 



Whiteoak & Green ( 1996 1 



3.2 Spectral analysis 

We extracted the spectrum for each XIS detector for the 
outer elliptical region, inner elliptical region and inner cir- 
cular region. The background has been extracted from cir- 
cular region, taking care not to include the 55 Fe calibration 
sources in the corners of the CCD. In Table 1, we gave de- 
tails of ellipses, circles and background (centre coordinates, 
sizes and angles). The spectra are grouped with a minimum 
of 50 counts bin -1 for region 1, 30 counts bin -1 for region 
2 and 25 counts bin -1 for region 3, and the \ 2 statistics is 
used. 

For spectral fitting of region 1, we first applied single 
temperature non-equilibrium ionization (NEI) model (VNEI 
in XSPEC; Borkowski, Lyerly & Reynolds (2001)) modified 



by Galactic absorption via the WABS multiplicative model 
( Morrison fc McCammon]|1983 1 . The parameters of the ab- 

electron temperature (kT c ) 
n e t) were set free while all 



sorbing column density (Nh), 
and ionization parameter (r = 
elemental abundances were fixed at their solar values (An- 



ders & Grevesse 19891. The model gave a large ionization 
parameter of r ~ 10 13 cm _3 s, which means that the plasma 
of the remnant has reached ionization equilibrium with a 
reduced \ 2 of 1.43 (948/664 d.o.f.). Therefore, we tried 
an absorbed collisional ionization equilibrium (CIE) plasma 
model (VMEKAL in xspec; Mewe, Gronenschild & van den 
Oord|(|1985|>;|Mewe, Lemen fc van den Oord|(|1986|);|Liedahl, 



Osterheld & Goldstein (19951)). In this case, the reduced \ A 
value (926/665=1.39) improved a little bit but it was clear 
that the model fitting needed a second component. There- 
fore, we added a second thermal component to VMEKAL, 
a CIE model (VMEKAL) and a NEI model (VNEI), re- 
spectively. From the model fitting of VMEKAL+VMEKAL 
we did not see two different temperatures, but when we 
thawed Mg, Si, S we obtained two different temperatures 
of ~0.6 and ~1.1 keV with a reasonably good reduced x 2 of 
1.15 (761.7/660 d.o.f.). However, we obtained unphysically 
low abundance of S and very high errors of abundances of 
Mg, Si and S. From the model fitting of VMEKAL+VNEI, 
we obtained two different temperatures of ~0.6 and ~1.4 
keV. When we thawed Mg, Si, S, the \ 2 value improved 
significantly (750.8/659=1.14). For both model fittings we 
used F-test to test whether the extra thermal component 
was a statistically significant addition. However, the F-test 
gave a probability of ~1.5 x 10~ 7 and ~6.9 x 10~ 9 that is 
the improvement is due to chance for VMEKAL+VMEKAL 
and VMEKAL+VNEI models, respectively. Then, we added 
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a power-law component to VMEKAL model. In this case, 
we obtained a reasonably good physical parameter values 
and their errors with an acceptable x 2 value (799/663=1.2). 
Then, we let Mg, Si and S vary which yielded a rather 
good reduced \ 2 value of 1.06 (700/660 d.o.f.). We also ap- 
plied F-test to our model fitting which gives a probability 
of ~1.2 x 10 -19 that the improvement is not due to chance. 
All these steps were repeated for region 2 and region 3. For 
region 2, VMEKAL model required a power-law component. 
However, region 3 did not require a second component. The 
XISO, XIS1 and XIS3 spectra of three regions were fitted 
simultaneously in 0.3—10 keV energy band and are given in 
Fig. 3 (a, b and c). The right panels of this figure show the 
individual components of the applied model, VMEKAL and 
power-law, for XIS1 only. Best-fitting parameters of applied 
models for three regions are listed in Table 2. All the errors 
are at 90 per cent confidence level. 

We also studied spectral variation across the remnant 
and gave the results in Fig. 4. To study the radial varia- 
tion of kT e and surface brightness we chose five elliptical 
annuli with 0.6 x 1.4, 1.3 x 2.1, 1.9 x 2.8, 2.7 x 3.6 and 
3.6x4.7 arcmin 2 regions centered at RA(2000) = 13 h 05 m 48 s , 
Dec. (2000) = -62°41'46" from XISO FOV. The spectra ex- 
tracted from the elliptical annuli were grouped with a mini- 
mum of 20 counts bin -1 . We applied our best-fitting model 
(an absorbed VMEKAL and power-law) starting from the 
central region within 0.6 x 1.4 arcmin 2 , we let absorbing col- 
umn density, electron temperature, normalization and abun- 
dances of Mg, Si and S vary. For all outer annuli, we applied 
same model fixing Nn to the value obtained from central 
region. 



4 DISCUSSION AND CONCLUSIONS 

In this work, we have studied X-ray morphology and spectral 
properties of G304. 6+0.1 using XIS data obtained with deep 
observation of Suzaku. As seen from Fig. 2, G304. 6+0.1 has 
a centrally filled X-ray morphology dominated by thermal 
X-ray emission, the radio contours of remnant are concen- 
trated at its shell. Morphologically, this remnant resembles 
MM SNRs. With the help of high quality spectral data we 
studied spectral characteristics of G304. 6+0.1 and compared 
our findings with the theoretical models proposed for MM 
SNRs below. 

The X-ray emission from the remnant is well described 
by a model of a thermal plasma that has reached ionization 
equilibrium and a non-thermal component (power-law) re- 
quired to describe the plasma state has a photon index (r) of 
~1.4. Dominant emission from this remnant is thermal (~90 
per cent of the total X-ray flux) and is described by the CIE 
plasma mod el (in CIE p lasma, n c t should be larger than 
10 12 cm~ 3 s ( Masai||l984 l) with an electron temperature of 
kT c ~ 0.75 keV. We have performed detailed spectral fitting 
for three regions. As seen from Table 2, spectral parameters 
for three regions are almost similar. We obtained a high ab- 
sorbing column density JVh in the range of (3.2 — 4.2) x 10 22 



cm which is consistent with the previous work of Combi et 



culations we used d=10 kpc. From the emission measure 
(EM = n c nnV , where n e is the electron density, rm is the 
hydrogen density, and V is the volume of the X-ray emit- 
ting plasma) we estimate the electron density of the plasma 
taking into account n = 1.2riH for the mean charge state 
with normal composition. Considering the possibility that 
less than the entire volume is filled, we calculate the volume 
of the X-ray emitting region from V — |7ri? 3 /, where R is 
the SNR radius and / is the filling factor. Then we carry the 
/ factor through our calculations to see the dependence of 
each derived quantity (for example the electron density, age) 
on /. We obtain the volume to be ~1.9 x 10 59 / cm 3 assum- 
ing the spherical radius of 4 arcmin at the distance d=10 
kpc and then the electron density to be ~2.3/ -1 ' 2 cm -3 . In 
calculating the age of the remnant we use the relation t/ti c , 
where r is the ionization age of the plasma. CIE model gives 
only a lower limit (10 12 cm~ 3 s) to this parameter and hence 
leads to a lower limit on the age t of the remnant which 
is calculated to be ~1.4 x 10 4 / 1/2 yr for G304.6+0.1. The 
mass of the X-ray emitting gas is calculated to be ^SSOf 1 ' 2 
M0 from M x = mun e V, where mn is mass of a hydrogen 
atom. As seen from Table 2, near/lower solar abundance val- 
ues obtained for each region and large X-ray emitting mass 
indicate that the X-ray emission results from swept-up am- 
bient medium as is generally the case in middle-aged SNRs. 
Relatively large n c (~2.3/ -1 ' 2 cm -3 ) and small size (~8 
arcmin) suggest that the remnant expands in a region with 
a large ambient density. 

For the non-thermal component that constituting ~10 
per cent of the total flux, our best-fitting photon index is 
in the range of 1.4—1.6, as seen in Table 2, which is consis- 
tent with that of classical young pulsars ranging in between 
1.1 and 1.7 ( Chakrabarty et al.|2001 |. But there is no com- 
pact or extended X-ray source reported at the location or 
in the vicinity of this remnant. Another possibility would be 
that the power-law emission is resulted from the synchrotron 
emission from shock-accelerated relativistic electrons, also 
as the Fermi detection of G304.6+0.1 ( |Wu et al.||2011[ ) ap- 
pears to be related to the remnant and not to a central 
compact object. However, our spectral fitting of each ellip- 
tical annulus gives different photon index value (r ~ 1.5, 
1.2, 1.3, 1.5), moreover the photon index obtained from the 
centre is unphysically large (F ~ 7.3) indicating that the 
power-law emission can not be related to the remnant itself. 
That our best-fitting model required no non-thermal com- 
ponent for region 3 may support this idea. Considering all 
these cases, it is still difficult to explain the origin of the 



non-thermal emission with Suzaku data. Combi et al. ( 2010 1 



also obtained non-thermal emission from this remnant with 
a photon index value in the range of 1.8—3.1 with XMM- 
Newton observation, which is softer than our photon index 
value, indicating particle acceleration in shock- fronts. 

In addition to morphological evidences, G304. 6+0.1 is 
a middle-aged (~10 4 yr) SNR, the plasma of the remnant 
is in CIE condition and it is located in an inhomogeneous 
medium. The general characteristics of MM SNRs have been 
given in Rho & Petre ( 1998 1 as (i) Members of this class have 



al. (20101, in the range of (2.5 - 3.5) xlO 22 cm~ 2 , that sup- 



ported the large distance d ^ 9.7 kpc adopted by |Caswell| 
et al.| ( |1975[) to the r emnant. Our Nn value is larger than 
that of Combi et al. (20101 therefore, throughout our cal- 



a centrally peaked morphology in the X-ray band character- 
ized by a thermal spectrum while a shell morphology in ra- 
dio, (ii) Their ages are expected to be ~10 4 yr, which means 
that the plasma should be in CIE condition, (iii) They are 
interacting with molecular clouds and/or HI clouds, which 
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means that they are located in a very dense medium. Consid- 
ering these properties we may say that G304. 6+0.1 is likely 
to be a member of MM SNRs as suggested by |Combi et al.| 
( 2010 1. The two basic scenarios explaining this class are; one 



is that the interior X-ray emission arises from the gas evap- 
orated from shocked clouds ( Cowie fc Mckee|1977 White & 



Long|[l99T I . The second one is that as an SNR evolves, the 
temperature and density of the hot interior plasma gradu- 
ally become uniform through thermal conduction and the 
temperature in the outer shells becomes lower and absorbed 
by the interstellar medium, the only detectable X-ray emis- 
sion emanates from the interior of the SNR (|Cox et al. 1999 



Shelton et al.|1999||Shelton, Kuntz fe Petre|2004| . To corn" 
pare our results with these models, we have studied radial 
variations of the kT B and the surface brightness as seen from 
Fig. 4. We have found that the temperature variation (kT c 
~ 0.66—0.81 keV) is almost uniform across the remnant as 
predicted by both models. On the contrary, surface bright- 
ness peaks at the center and declines towards shell regions 

0.03X10" 12 
in 0.3—10 keV, which is consistent 
with evaporating clouds model. G304. 6+0.1 is interacting 
with an OH maser source at 1720 MHz ( |Frail et al.]|1996[ |. 
Recently Lee et al. (20111 suggested that the remnant was 



from ~2. 34x10 12 erg s 'cm 2 arcmin 2 to 
erg s _1 cm _2 arcmin -2 



interacting with nearby molecular cloud. Furthermore, |Wu 



et al. (20111 have detected GeV emission from the SNR 



and a number of 7-ray sources in its vicinity. The thermal 
conduction model requires a relatively high density ambient 
medium, while the evaporation model requires dense clouds. 
In this respect, G304. 6+0.1 is consistent with both models. 
We present here the results from the Suzaku archival 
data of G304. 6+0.1. The deep XIS spectra consist of two 
components; thermal emission coming from the plasma is in 
CIE condition with kT c ~ 0.75 keV, a high Nu value and 
the non-thermal emission has a photon index of ~1.4. This 
value is consistent with that of classical young pulsars, but 
there is no observed compact source reported. So, the origin 
of the non-thermal emission is still unclear. Thermal emis- 
sion originates from the shocked ambient medium with an 
electron density of n c ~ 2.3/ -1 ' 2 cm -3 . The morphological 
and X-ray spectral characteristics imply that G304. 6+0.1 is 
a member of MM class. 
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Table 1. Centre coordinates, sizes and angles of region 1, region 2, region 3 and background. 



Regions Centre Coordinates Size Angle 

RA(2000) , Dec. (2000) (arcmin) (degrees) 

/ h m s o / //\ 



1 

2 

3 

Background 



13 05 44 
13 05 55 
13 05 35 
13 04 52 



-62 41 54 
-62 41 11 
-62 43 33 
-62 37 16 



3.2 x 4.1 
1.9 x 3.4 

1.24 

1.56 



354 
342 







Table 2. Best-fitting parameters with corresponding errors at 90 per cent confidence level in 0.3- 
energy band. 



10 keV 



Component 



Parameters 



Region 1 



Region 2 



Region 3 



Wabs 

VMEKAL 

Abundancq a 



Power-law 



AT H (xl0 22 cm- 2 ) 

fcT c (keV) 

Mg 

Si 

S 

norm (photons cm" 2 s 

eaM 

FIujF] 

Photon Index 
norm (xl0~ 2 photons cm 
Fhrsj^ 

X 2 /d.o 



- 1 ) 



- 1 ) 



3.9±0.1 
0.75 ±0.01 

1.3 ±0.2 
0.9 ±0.1 
0.6 ±0.1 

0.93 ±0.06 
11.1 ±0.8 
4.6 ±0.1 

1.4 ±0.1 
1.17 ±0.23 

5.5 ±0.2 
700/660=1.06 



4.2±0.1 

0.74 ±0.01 

1.2 ±0.2 

0.9 ±0.1 

0.5 ±0.1 

0.65 ±0.11 

7.7 ±1.1 

3.2 ±0.2 

1.6 ±0.2 

0.46 ±0.21 

3.6 ±0.2 

747/657=1.14 



3.2±0.1 
0.84 ±0.03 
0.9 ±0.2 
1.1 ±0.1 
0.4 ±0.1 
0.07 ±0.01 
0.8 ±0.1 
0.6 ±0.1 



253.8/192=1.32 



a Abundances are relative to the solar values ( Anders & Grevesse|l989 1 . 

Emission measure EM = I n c nnd~V in the unit of 10 5y cm -3 , where n e and nn are number densities of 
electrons and protons, respectively and V is the X-ray emitting volume. 
c Absorption-corrected thermal flux in 0.3 — 10 keV energy band in the unit of 10 — lx erg s"'cm -2 . 

Absorption-corrected total flux (VMEKAL plus power-law components) in 0.3 — 10 keV energy band in 
the unit of 10 -11 erg s _1 cm -2 . 
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Figure 1. Suzaku XISO image of G304. 6+0.1 in 0.3—10 keV energy band. The outer and inner ellipses represent region 1 and region 
2 respectively and inner circle represents region 3. The dashed circle and dashed square show background region and FOV of XISO, 
respectively. 
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Figu re 2. Suzaku XISO im age of G304. 6+0.1 in 0.3—10 kcV full energy band overlaid with contours of the 843 MHz radio image from 
|Whiteoak &: Green| |l996^ . The X-ray image is smoothed with a Gaussian kernel of <r=2 arcsec, and units are in counts s — 1 . Overlaid 
contours arc spaced linearly in intensity of 0.07, 0.26, 0.44 and 0.62 count pixel -1 . 
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Figure 3. Upper left panels: Suzaku XIS (XISO:black, XISl:red and XIS3:grccn) spectra have been fitted by an absorbed 
VMEKAL+power-law models in 0.3—10 keV energy band simultaneously, a) region 1 and (b) region 2. Upper right panels: The red 
curve represents the best-fitting of two component model (VMEKAL+power-law). The individual components are represented by the 
cyan curve for the VMEKAL model and the yellow curve for the power-law model for only XIS1 for region 1 and 2. Bottom panel (c): 
XIS spectra have been fitted by an absorbed VMEKAL model in 0.3—10 keV energy band simultaneously for region 3. The lower panels 
show the residual in units of a. 
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Figure 4. The radial profiles of electron temperature (left panel) and surface brightness (right panel) obtained from XISO spectra. 



